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Abstract 

We study the effects of Lorentz invariance on relativistic constituent quark model wave functions. The 
model nucleon wave function of Gross et al. is constructed such that there is no orbital angular momentum 
and that the spin-dependent density is spherical. This model wave function is claimed to be manifestly 
covariant. We consider two possible interpretations of the nucleon wave function in an arbitrary reference 
frame. In the first, the seeming covariance of the matrix elements of the electromagnetic (em) current 
arises from using the Breit frame. Matrix elements have a different appearance in any other frame. In 
the second interpretation, the em current is covariant yet it is not consistent with the general structure 
required by QFT, e.g. the wave function of the incoming nucleon depends on the momentum of the 
outgoing nucleon and vice versa. 
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I. INTRODUCTION 



Recent Jefferson Laboratory data on em form factors of the nucleon have created much theoret- 
ical interest. The key finding is that the ratio of the proton's Ge/Gm falls rapidly with increasing 
Q 2 [l, 2]. But new results for the neutron electric and magnetic form factors have been or are about 
to be obtained; see the reviews [3-5]. 

It was argued[6] that reproducing the measured ratio Ge/Gm ratio requires a relativistic treat- 
ment that includes the effects of the quark's non-zero orbital angular momentum. Ref. [7] intro- 
duced the idea of using the rest-frame to rest-frame matrix elements of a spin-dependent charge 
density operator to exhibit the influence of the orbital angular momentum. In particular, (for a 
model without explicit gluons) the probability for a quark to have a given momentum, K, and a 
given direction of spin, n is given by [7, 8] 

/<7 3 r - 1 

7^ elK >(r)0-(7° + 7 -n 7 5M0), (1) 

where O is Q/e, the quark charge operator in units of the proton charge for the spin-dependent 
charge density or O = 1 for the spin-dependent matter density. The matrix element of the 
operator pe>(K, n) gives the spin-dependent matter densities. The quark field operators ^(r), i/>(0) 
are evaluated at equal time. The rest-frame matrix element of this density operator in a nucleon 
state of definite total angular momentum defined by the unit vector s, \^ s ) is 

p (K,n,s) = (* B |p (K,n)|* 8 ), (2) 

where the subscript O = Q, or O = 1 specifies the operator used in Eq. (1). The most general 
shape of the proton, obtained if parity and rotational invariance are upheld can be written as 

Po (K, n, s) = A (K 2 ) + B (K 2 )n • s + C (K 2 ) (n • K s • K - • K K 2 ^) , (3) 

with the last term generating the non-spherical shape. The effects of non-vanishing orbital an- 
gular momentum cause the matrix elements of the spin-dependent density operator Eq. (2) to 
be non-spherical and yield a non-zero value of the coefficient Ce>(K 2 ). While no experiment has 
been constructed to measure the spin-dependent density, this quantity can be evaluated using the 
techniques of Lattice QCD, and has been measured in condensed matter physics[9]. 

Gross & AgbakpeflO] constructed a relativistic constituent quark model that was claimed to have 
a spherical shape. However, these authors did not consider the spin- dependent density operator. 
When we [8] used the wave function of [10] to evaluate the matrix element of the spin-dependent 
charge and matter density operators, a non-spherical nucleon shape was obtained. More recently, 
Gross et al. [11] claimed to find a covariant constituent quark-diquark model that describes all 
the available em form factors, but has no orbital angular momentum. In this case, the shape of 
the proton as determined by the rest-frame matrix element of the spin-dependent density matrix 
is indeed spherical. The question of whether or not it is possible to find a covariant model that is 
pure S-wave is an interesting one that we examine here. 

We next describe the wave function of [11], using the notation of that reference. The nucleon 
wave function ^N(P,k), of total four-momentum P and di-quark four-momentum k, is given by 
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the expression 
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]=MP,k) 0} 75 f P u(P,s), 



V2 



(4) 



which is a sum of contributions from a spin-isospin (0,0) diquark and a spin-isospin (1,1) diquark 
and ^0,1 are Lorentz scalar functions. The polarization vectors ep are given by the expression 



where is a genuine relativistic polarization vector of a vector particle (di-quark in the present 
case) 6^ • k — 0. This quantity is denoted by r\ = in [11]. The operator Op is a Lorentz 
transformation, with 



the four-momentum state (E s , 0, 0, k) to the diquark rest frame (m s , 0, 0, 0), and finally Bp boosts 
the vector (M, 0, 0, 0) to the moving frame (E P , 0, 0, P). The wave function \I> N satisfies the Dirac 
equation because P commutes with 7^p. As stressed in ref. [11], the essential difference between 
this model and the one introduced in Ref. [10] is that in the nucleon rest frame, the wave function 
(4) contains absolutely no angular dependence of any kind. 

We discuss the general requirements for covariance and a proper treatment of a relativistic 
constituent (quark di-quark) model. We study two interpretations of [11] based on two different 
generalizations of the boost Bp in Eq. (4) to the case of the arbitrary P and find that, using the first 
(conventional) interpretation, the model wave function of [11] does not satisfy these requirements 
because it is not covariant, and as a result produces an em form factor which is not Lorentz 
invariant. Using the second (unconventional) interpretation leads indeed to a Lorentz invariant 
em form factor. However, the use of the second interpretation does not yield a model that satisfies 
the basic requirements (that any composite quark model must satisfy) discussed in Sect. IV. We 
summarize in Sect. V. 

II. COVARIANT VECTOR DI-QUARK WAVE FUNCTION 

Let us denote the vector di-quark wave function as typ jS (k,e), defined as 



where (P, s\ and \k, e) are nucleon and di-quark eigenstates and q(0) represents a quantized quark 
field operator. Note that the dependence on the polarization vector e and nucleon spin is made 
explicit. Lorentz invariance requires that 



Ep — Op€i~, 



(5) 





(7) 



<j> P>s (k,e)~U(P, S )r^P,k)e 



k- 



(8) 



where r M (P, k) is a covariant vector : 



r M (p, k) = A lhllx + p 75 /c M + c 75 p m + • • • , 



(9) 
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where A, B, and C are Lorentz scalar functions built from the four vectors P and k. The forms 
Eq. (8), and Eq. (9) have been known for a long time [12] and have been applied recently [13, 14]. 

To see how this Lorentz invariance of Eq. (8) works in practice, consider the relevant particular 
example of the matrix element of the em current: [15] 

<P + ,s'\J a (0)\P^,s >=M a = I d 4 k J2*P + Ak,eh a Vp-,s(k,t) 

~ f d'k ■■■Y,U(P + ,s')T,(P + ,ky kl a elT u (P^k)U(P-,s), (10) 

where initial and final nucleon four- momentum are denoted as P_ and P + . In Ref. [11] (their 
eq. (11)) P± are explicitly chosen in the Breit frame. Here the only restriction is that P + = P_ + q 
and P 2 . = M 2 , where q is the four-momentum of the virtual photon and M is the nucleon mass. 
The quantity Ai a should be explicitly Lorentz invariant. The sum over polarizations is performed 

as 

E*! = ^— (id 

e 

where m is the di-quark mass. Thus one finds 

/k^k v 
d A k ■■■ u(p + , s')r M (P+, k)r(— - ^)r,(p_, k)u(p_, s). (12) 

The result Eq. (12) has a manifestly covariant form as a Lorentz four-vector that results from the 
use of Eq. (11). 



III. NON-COVARIANT WAVE FUNCTION OF [11] 

We need the nucleon wave function Eq. (4) in an arbitrary reference frame. For this we use 
first the straightforward (conventional) definition of the boost Cp which is explicitly stated in [11], 
between eqs (6) and (7): "The spin states are analogous of (2) and (3), and their form in an 
arbitrary frame is obtained by boosting the nucleon to momentum |P| = P along the z direction 
and then rotating". This is equivalent to the statement that the z-direction cannot be favored 
among others in Eq. (4). We will show that such model corresponding to the "first interpretation" 
of [11] is not consistent with Lorentz invariance. We note in advance that the essential point will 
be that different polarization vectors ep + and ep_ enter into the sum over polarization vectors. 

Here consider the following (first interpretation) generalization of the vector-diquark part of the 
nucleon wave function Eq. (4) to an arbitrary reference frame: 

* P (fc) ~ U(P, s) 7 57m4 = U(P, s) lblll {CpC- k l e k y 

C?e k = (p,?) = eo (13) 

where P and k are the nucleon and diquark on mass shell momenta, C k 1 is the boost transformation 
C-^k = (Vk?, 0), €k is a genuine relativistic polarization vector of a vector particle (diquark in 
the present case) e& • k — 0, which is denoted by r\ = e& in [11]. Furthermore, eo is the diquark 
polarization four-vector in the diquark rest frame. Our notation here differs slightly from that of 
[11] because (we use the notation C k instead of R k Bk) and because the quantity C k l is not exactly 
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the same as B^ l R~^ 1 of Eq. (6): for the sake of simplicity we not make the effects of the rotation 
explicit. This simplification does not affect our conclusions [16]. 

Let us emphasize that the first interpretation consists in replacing the boost Bp (with z directed 
P) in the Eq. (4) by the boost C P (with arbitrary P). As noted above in Eq. (9), Lorentz 
invariance requires r M (P, k) to be a covariant four- vector in any quark-diquark wave function 
'fyp(k) ~ U(P, s)r /i (P, k)e%.. The result Eq. (13) is not consistent with this requirement because 
the quantity 757^ (CpC^^ is not a four- vector. In particular, the explicit appearance of the 

product of boosts, Cp,C^ x , breaks covariance. Neither Cp nor C^ 1 is a covariant tensor. To see 
this we derive the boost tensor from the expression for the boost Eq. (2.8) of [17] (note that a sign 
misprint in that equation is fixed here): 

^ = * " WT^ k - k " ~ + W^ s -' k " ~ fT^°" < 14 > 

The result Eq. (14), with the explicit presence of the index 0, makes it clear that (jC^ 1 ) 1 ^ is not a 
covariant tensor. Similarly 

where we have used P M = 25qP° — P M . 

Lorentz invariance is lost if one uses the wave function of [11] ^p(k) ~ U(P, s)T ^{P,k)e P 
because the sum over diquark polarization, J2 € e P + e p-i that enters in the matrix element of the em 
current is not Lorentz invariant. Let us calculate the diquark polarization sum D^ v '. We find 

Dp, = $>P + yepJ„ = ^P + )> a (Cp_te , = {C P+ %{Cp_)l (16) 

e a/3 i=l,2,3 

Note that here the quantity D^ U (P + , P~) is seen to be a sum of product functions, with one 
function depending only on P + and the other depending only on P~ . 

We use the expression for the boost Eq. (15) to evaluate Eq. (16), with the result 

M 2 — P + ■ P_ 
" <VoO*) (p o + M) (p o + M) 9,u 

PyPy Py Pji P^Pv P+ ' P- 

' 1 x-.n 71 71 r ~r / 7-.n 71 r \ 71 r ' 
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(P° + M)M (P°+M)M (P° + M)(P° + M)M' 

^ Prf P+S^P+P- + PjM) + SjP-(P+P- + P°M) 
po + M P°+M (P° + M)(P° + M)M ' [ ' 

A brief inspection shows that D^ u , as obtained in a general reference frame, involves the non- 
covariant expressions 5 M o as we U as explicit three-vectors and therefore is not a covariant tensor. 
This result means that the wave function of [11] is not covariant and that the expressions for 
matrix elements of the em current that result from using Eq. (17) are not covariant. 

However, one can be fooled by using one particular frame- the Breit frame. In this case, the 
four-vectors P± are given by P+ = (E, 0, 0, Q/2), P + = (E, 0, 0, — Q/2). It is also useful to note 
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that the non-covariant expression 5q = (1, 0, 0, 0), can be written in an apparently covariant form 

fP-4- + P Y 
y/AM 2 - (P + - P_) 2 

We proceed by evaluating Eq. (17) in the Breit frame. Use Eq. (18) and P + = — P_ = P, 
P+ • P_ = 2P 2 - M 2 in Eq. (17) to obtain 

9 M 2 — ?P 2 

D =(p+ + p-)(p+ + p-) =12 a 

H u \ r » H I \ v ^ r v > ^pQ _|_ M) 2 Pq ^ 

, + + p;k (m 2 - p^ 



(P° + M)M (P° + M)M (P° + M) 2 M 2 



2P°(P° + M) 2P°(P° + M) 
[P+(P+ + P~) + (P+ + P-)P-][2P 2 - M 2 + P°M = (P° + M)(2P° - M)] 

2(P° + M) 2 P°M 



2P 2 M 2 

1 P+P" 

= (^ + + p;)(p? + ^) M2 + P+ . P _ - ^ - -^r ■ (19) 

This result, obtained previously in Ref. [11], has a illusory covariant appearance, resulting from 
the explicit use of the Breit frame. The expression Eq. (19) would not be correct in any frame other 
than the frame where 3D parts of P + and P_ are collinear. In particular, the factor M 2 + P + • P_ 
that appears in the denominator of Eq. (19) violates the sum of product functions form of Eq. (16). 



IV. WAVE FUNCTION WITH UNCONVENTIONAL POLARIZATION VECTORS 

We found out above that the first interpretation wave function of [11] is not covariant, and as a 
result it fails to produce a Lorentz invariant em form factor. Then we can ask the question, what the 
wave function of [11] should be in the arbitrary reference frame to produce the Lorentz invariant em 
form factor given that the z direction moving nucleon is described by the wave function of Eq. (4)? 
The unambiguous solution to this question is to use the "unconventional" diquark polarization 
vectors £(P±) 1 in the covariant wave function (8) and correspondingly in Eq. (10) instead of usual 

£(P + ) = Ae(Z + ) = AC z+ e , £ T (P„) = e T (Z„)A T = e T C T z _k T , (20) 
where A = C P++P ^ is the boost transformation defined as 

A-\P + + P_) = ( V / (P+ + P-) 2 , 0) = Z + + Z_, Z ± = A-'Pi. (21) 



We understand that these polarization vectors are suggested (and named "unconventional") by F. Gross [18] for 
the construction of the pure S-wave wave function. 
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— * — * 

So the 4- vectors Z + and Z_ are collinear, Z_ = —Z + . Indeed Lorentz invariance of the em current 
implies that it transforms under the Lorentz transformation L as 2 



< P + ,s'\J a (0)\P^,s >=M a = U(P + ,s')M a (P + ,P_)U(P_,s), 
M a (LP + , LP_) = L«S L M?(P + , P_)S+, where S Ll a S+ = {L~%y 



(22) 



Therefore the em current can be boosted to the arbitrary frame as follows 



M a (P + ,P_) = M a (AZ + ,AZJ) = A%S K M p (Z + ,Z_)St, 



(23) 



where 



M' 




(24) 



is expressed in terms of the wave functions, (4), of nucleons moving along z direction. Simple 
algebra leads unambiguously to the unconventional polarization vectors of Eq. (20): 



where Lorentz invariance for the scalar functions was used, ipi(Z±,k) = ipi(P±, Ak). Now, using 
(25) in the first line of the Eq. (22), the expression for the em form factor in terms of the new 
(unconventional) polarization vectors reads 



where T' U (P_, k) = 757^1 {P-i k) corresponds to the wave function (4). Below the primed T' U (P_, k) 
denote more general matrices used in the Eq. (32), i.e. in the wave functions of [11] corresponding 
to the round spin-dependent density. They should not be confused with the (unprimed) T V {P_, k) 
used in the covariant wave functions (8), which according Eq. (9) are strictly covariant and do not 
depend on the momentum of other nucleon to be consistent with the well known principles. 

To assess the expression (26) we need to compare it with Eq. (12), the most general expression 
for the nucleon current, < P + , s'\ J a \P-, s >, in the relativistic quark-spectator diquark model 
(with diquark on mass shell) with the photon interacting only with the quark (using a point-like 
photon-quark interaction and ignoring explicit factors of charge), where T^(P + ,k) and T u (P^,k) 



2 Just the same arguments of Lorentz invariance for the quark di-quark wave function (7) unambiguously lead to 
the covariant form (8). This paper addresses possible violations of general principles if forms inconsistent with 




(25) 




(26) 



(8) are used. 
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are covariant functions of Dirac 7-matrices and the momentum variables. The factor — % q' 11 ' 

arises from the axial-vector di-quark propagator and must present in all quark spectator axial- 
vector di-quark models . 

We stress that (12) must hold independently of the choice of the particular form the diquark 
polarization vectors. A given model is defined only by the choice of a specific form for r M (P + , k) 
and r„(P_, k), the rest is fixed by (12), there is nothing left to choose. The results of [11] or any 
other model must be consistent with Eq. (12), they should correspond to particular choices of a 
covariant F V (P±, k). 

It is useful to discuss the result (12) in the operator formalism of the quantum field theory. We 
compute the matrix element < P + |g7 Q g|P_ > by inserting the sum over a complete set of states: 

M a =< P+, s'\q"f a q\P., s >~ < P +' s '\v\ n >l a < n\q\P-, s > . (27) 

n 

Truncating the sum to the term of the valence quark model, \n >< n\ — » f^2\k,i >< k,i\, 
(where \k,i> are diquark states) gives: 

M a ~ / al 4 k5 + (k 2 - m 2 ) U(P+, s')f M (P+, k)-f a ^(k, i)e v {k, i)r„(P_, k)U(P-, s) (28) 
J i 

where the diquark "conventional" polarization vectors, e M (A;, i), e u (k, i), are exposed. The integrand 
of Eq. (28) is a sum of contributions factorized in the nucleon momenta, P + ,P_, it is a sum of 
terms each of which is of the form of product of scalar functions:Pi(P + , k)F2(P-, k). This property 
that the integrand is a sum of products is a signature property of a valence quark model (covariant 
and non- covariant). 

Thus the relevant requirements for obtaining a correct evaluation of a relativistic valence quark 
spectator axial di-quark model are summarized in Eq. (12) and Eq. (28). 

We find that D IXV {P + ,P_) in Eq. (26) is indeed the covariant tensor given in (28) of [11]: 

D»"(P + ,P-) = £«i>+) x e(P-) = Ae(Z + ) x e T (Z_)A T 

= a {(z; + z- )( zt + z-) M2 + \ + z _ -^- z -i§~Y 

1 P+P- 

= (p+ + p;) (p+ + P ») M 2 + P+ . P _ - 9r - -5^ ( 29 ) 

But a brief inspection shows that the result (26) is not consistent with the general form (12). The 
integrand of the general requirement (12) is of the form of a sum factorized of terms of the form 
Pi(P + , k)F2(P-, k) in the nucleon momenta, P + ,P_, the integrand of the expression Eq. (26) is 
not factorisable due to the presence of the denominator M 2 + PP_ in D flu (P + , P_). There is no 
way to derive the integrand of (26) from (12). 

We explain this in more detail by explicitly using the model of [11] in Eq. (26). This is to 
illustrate that the Dirac operators do not yield a factor of M 2 + P + ■ P_ in the numerator that 
cancels the one in the denominator. The model takes the form r'(P_,fc) = 7„7 5 -?/>i(P_, k), with 
■01 a scalar wave function. Use these T^s and Eq. (19) in Eq. (26) to obtain: 

M a = J d 4 k5 + (k 2 -m 2 )^ 1 (P + ,A;)^i(P_ ! A;)f/(P + ,/)7^(P + ,P_)7 Q 7^(P_, S ) 
= J d A k5 + {k 2 -m 2 )4> 1 (P + ,k)MP-,k)U(P + ,s f ) 



4M (P» + Pf ) a 
M 2 + P+ • P_ 7 



U(P.,s). (30) 
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The explicit appearance of a term inversely proportional to M 2 + P + ■ P_ shows that the model [11] 
violates the general requirements expressed in Sects II and IV. The use of the covariant expression 
(19) for D^ u leads to a contradiction with well-known principles. 

It is important to see that the transformation A and vectors Z + and Z_ depend on both 
momenta, P + and P_ and thereby so do the unconventional polarization vectors (therefore the 
notation £(P±) with only one momentum argument is misleading)! This fact is not emphasized in 
[11]. It might be interesting to show that this dependence on two momenta is completely contained 
in Wigner rotations, Wj(P±, Z±), usually involved in the Lorentz transformation of polarization 
vectors 

£(P+,i) = ke\Z + ) = Y^e j (AZ + )W;(AZ + ,Z + ) = ^(P + )Wj(P + ,Z + ) (31) 

3 3 

where the polarization vector, e J (P+) = Cp + e J Q , depends only on P + just in the same way as the 
genuine polarization vector e 3 k from (8) depends on k. The whole dependence on the two momenta, 
P + and P_, is contained in the Wigner (purely 3D) rotations, Wj(P+, Z + ), through Z + = A _1 P + . 
As a result the wave function of the initial (final) nucleon in the Eq. (26), 

* P _~Z"(P_,{)>fr' v (P-,k)U(P-), * P+ ~ U{P + )T'^P + ,k)e{P + ,i), (32) 

depends not only on the momentum of the initial P_ (final P + ) nucleon, it depends on both, P + 
and P_, which is beyond the very idea of quark models. This dependence on the both momenta 
comes via the "unconventional" diquark polarization vectors, they depend on both the initial and 
the final nucleon momentum, P + ,P_. 

Having shown that the wave function of [11] corresponding to Eq. (4) is not consistent with 
well-known principles, here we can discuss restrictions these principles impose on the general form 
of T^(P + , k) if one is to use the wave functions (32). 

The minimal consistency condition is 

E ^( P +> s ')f M (P+, k)^(k, i)<*(k, i)r„(P_, k)U(P-, s) 

i 

= J2 U(P + s')T^P + , k)e(P + , i)C(P-, ih a K(P-, k)U(P_s) (33) 

which restricts T' v because r^(P_, k) is restricted to being covariant and independent of P+. This 
leads to 

k^k u 

u{p+, s')r^(p + , k^i—r - -rWP., k)u(p_, s) 

= U(P + s')T^P + , k)D^(P + , P-)rK(P-, k)U(P.s). (34) 

Although Eq. (34) admits a covariant solution for T' u (P-,k) it requires T' u (P-,k) to depend also 
on P + . A solution to Eq. (34) is 

rt(P-, k) = (A^A-^r^p., k), ryp + , k) = f Q (p + , k){C k RC- z \\-X, (35) 

where R is an arbitrary 3D rotation. The mentioned dependence on the momenta of the both 
nucleons comes from CkRC^-hr 1 . The solution (35) can be verified by using the identity 

^4- - 9 aP = {C k RC- z \K-XD^{P + , P-){C h R£?_K-y v (36) 
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The arbitrariness of R follows from the identity 

The expressions in Eq. (34) are closely related to the generalized parton distributions (GPD) 
corresponding to the parton (quark) momentum P_ — k. We would like to note that for any given 
model of the nucleon one can construct such GPD related function[19] which reproduces exactly 
the em current of the given model and at the same time it corresponds to the round spin-dependent 
matter density. The problem is that in such GPD the initial and final nucleon momenta are not 
factorized and therefore it cannot be explained in the framework of a valence quark model. 

V. ASSESSMENT 

Using the first of two interpretations discussed above (conventional polarization vectors), we 
have shown that the seemingly covariant appearance of the expressions of [11] results from the 
explicit use of the Breit frame. This failure to maintain covariance results from using the polariza- 
tion vectors ep instead of to describe the vector di-quark wave function. However, this is a very 
important point in the present context because it is exactly the use of ep that allows the construc- 
tion of a model wave function without orbital angular momentum. As noted in Ref. [11], the result 
(19) has no angular dependence, so the evaluation of the matrix element of the spin-dependent 
density operator would yield a spherical shape. However, this roundness is caused solely by the 
lack of Lorentz invariance. Using the polarization vector would lead to a model much like that 
of [10], which does have a non-spherical shape, as measured by the spin-dependent matter density. 

Given the importance of Lorentz invariance, we have derived the general form of the wave func- 
tion that would produce a Lorentz invariant em form factor. This derivation led us unambiguously 
to the "unconventional" diquark polarization vectors £(-P±), suggested in [18], to interpret the 
wave functions used in [11]. These yield covariant results. Unfortunately neither the derived wave 
functions nor the expression for the em current the satisfy well-known principles discussed in the 
Sect. IV. For example, the wave function of the incoming nucleon depends on the momentum of 
the outgoing nucleon too and vice versa, so the em current cannot be written in the traditional 
form of the convolution of two proper wave functions. 
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